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We present a calculation of the wavevector-dependent subband level splitting from spin-orbit 
coupling in Si/SiGe quantum wells. We first use the effective-mass approach, where the splittings 
are parameterized by separating contributions from the Rashba and Dresselhaus terms. We then 
determine the parameters by fitting tight-binding numerical results obtained using the quantitative 
nanoelectronic modeling tool, NEMO-3D. We describe the relevant parameters as a function of 
applied electric field and well width in our numerical simulations. For a silicon membrane, we find 
the bulk Rashba parameter to be linear in field, a — a^E z with a 1 ~ 2xl0~ 5 nm~ 2 . The dominant 
contribution to the spin-orbit splitting is from Dresselhaus-type terms, and the magnitude for a 
typical flat SiGe/Si/SiGe quantum well can be as high as 1/xeV. 
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Silicon is a leading candidate material for spin-based 
quantum information processing^* 2 .. Its spin-orbit cou- 
pling (SOC) is relatively weak and the hyperfine cou- 
pling can be eliminated by isotopic purification-. This 
means that the spin lifetimes should be long. One 
way to measure a spin lifetime is to use electron 
spin resonance (ESR) on two-dimensional electron gases 
(2DEGs) in silicon quantum wells (QWs)*i2j&, where 
the D'yakonov-Perel— mechanism accounts for the re- 
laxation. SOC may also be measured directly, using 
EST— £ or photocurrentsiS. In order to compare with 
experiment, however, we need the wavevector-dependent 
SOC Hamiltonian, which must be calculated atomisti- 
cally. That is the aim of this paper. We shall focus 
on Si layers grown in the [001] direction with Si^Gei-^ 
layers on either side. 

Most calculations of the DP relaxation have used the 
Rashba Hamiltonianii, which is of the form 



Hr 



Oy ky) 



(2) 



H R = a (E X ,N, v) (a x k y - a y k x ) . 



(1) 



<Tj are the Pauli matrices and fcj are the two-dimensional 
wavevector components. N is the number of atomic Si 
layers in the well. v is the valley degree of freedom. 
We focus on the lowest electric subband, in which case 
the valley degree of freedom is two-valued- 1,2 and v is a 
two-by-two matrix. 

a(E z ) depends (in lowest order) linearly on E z , the 
external electric field, and often only this term is kept. 
E z is of order 1 — 5 x 10 7 V/m in heterostructures or 
MODFET devices^. The large magnitude of the field 
makes it important to examine the linearity assumption, 
and that is one of the purposes of this paper. The Rashba 
is generally thought of as a bulk effect. However, de 
Andrade e Silva et al. pointed out that surface effects 
may also be important, and this assumption should also 
be re-evaluated. In addition, when a detailed treatment 
of the surface effects was done by Nestoklon et al. in the 
absence of an applied field E z , they showed that in this 
case one obtains a term 



if N is od d 15 ' 16 . This is a Dresselhaus-like termi 7 - in 
that it arises from inversion asymmetry. The lowest- 
order term in E z for (3 (E z , N,v) is a constant, so that 
if N is odd then we get a surface-induced spin splitting 
even in the absence of an external field. 

The symmetry considerations refer to an ideal free- 
standing Si layer, or a layer that is sandwiched between 
two identical Si^Gei^a; layers that are treated in the vir- 
tual crystal approximation. Real Si^Gei-^ layers have 
substitutional disorder that destroys all symmetries. 

The purpose of this paper is to determine the func- 
tions a (E Z ,N, v) and (3 (E z , N, v) both for the ideal, free- 
standing case and for the case of a well confined in a 
Si/SicGei-a; heterostructure. These functions determine 
the spin properties of electrons in Si quantum wells. We 
shall focus particularly on the question of which term 
dominates for quantum wells with realistic values of E z 
and N. 

The numerical tight-binding calculations are per- 
formed with NEMO-3D 18 on nanoHUB.org computa- 
tional resources^. In NEMO-3D, atoms are represented 
explicitly in the sp 3 d? s* tight-binding model, and the va- 
lence force field (VFF) method is employed to minimize 
strain 20 . NEMO-3D enables the calculation of localized 
states on a QW and their in-plane dispersion relation 
with a very high degree of precision, allowing to extract 
the splittings along the in-plane directions in k space. 
Note that earlier calculations of a and j3 have used the 
virtual-crystal approximation for the interfaces. Since 
both of these quantities are dominated by atomic-scale 
interface effects, this is a rather crude approach. 

This paper is organized as follows: Sec. [IT] discuss the 
symmetry operations of a silicon membrane. In Sec. IIII1 
expressions for the SOC in a ^-functional effective mass 
approach are given. We present a qualitative picture in 
Sec. IIV1 Sec. IVl contains the numerical results for ideal 
Si QWs and Sec. ED the ones for a SiGe/Si/SiGe het- 
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erostructure. We conclude in Sec. I VIII with a summary 
of the results obtained. 



II. SYMMETRY 

We give here for clarity the symmetry operations of 
this system in the ideal case, since they are far from ob- 
vious. We stress that this is only an expanded discus- 
sion of the analysis already given by Nestoklon eiaL— . 
The lattice considered as a bulk sample has the dia- 
mond structure with a tetragonal distortion due to the 
Si^Gei-a; layers: the [001] axis along the growth di- 
rection is compressed relative to the in-plane [100] and 
[010] axes. In a symmorphic lattice this would simply 
reduce the point group symmetry from the cubic group 
Oh, with 48 operations, to the tetragonal group D^h, with 
16. (Recall that a symmorphic space group is one that 
is generated by translations and by rotations and reflec- 
tions about a point. A nonsymmorphic lattice requires 
combined operations such as screw axis and glide plane 
operations in its generating set.) For the diamond lat- 
tice, a nonsymmorphic lattice, there is no point group, 
but the factor group still has 16 operations, realized as 
follows. There are 4 proper rotations through 180° using 
the 2-fold [100], [010], and [001] axes passing through the 
origin, which may be taken at the position of any atom. 
In addition there are two ±90° improper rotations about 
the [001] axis: the rotation is followed by a reflection in 
the (001) plane. There are also reflections with respect 
to the (110) and (110) planes. Any of these 8 opera- 
tions may be combined with an inversion through a point 
midway between any pair of nearest neighbors. When 
we consider a layer, all 16 operations might appear to 
preserve the positions of the interfaces, whose presence 
would therefore not reduce the symmetry. However, this 
turns out to be not quite the case. 

For a layer with an odd number N of atomic lay- 
ers, take the origin (0, 0, 0) at an atom in the central 
plane. Then the other atoms in this plane lie at the 
points (a x /2) [m (1, 1, 0) + n 2 (1, —1, 0)], where a x is the 
in-plane lattice constant and n,\ and n 2 denote any in- 
teger. Each atom in the central plane is accompanied 
by another one in the plane z = a z /4 shifted from it 
by (a x /4, a x /4, a z /A) , where a z is the lattice constant in 
the z-dircction. The atoms in the plane at z = — a z /4 
are shifted with respect to the atoms in the central plane 
by (ax/4, —a x /4, — a z /4) . The atoms in the planes at 
z = ±a z /2 are shifted with respect to the atoms in 
the central plane by (a x /2, 0, ±a 2 /2) . The positions of 
other atoms can be found by translating this layer of 
thickness a z by integer multiples of (0,0, a z ) , so to un- 
derstand the symmetry properties of the full layer it is 
sufficient to consider a layer with N = 5. Simple re- 
flection in the z = plane is not a symmetry, since it 
interchanges the z — a z /4 and z = —a z /A layers, whose 
in-plane shift is (a x /2) (0, 1,0) . The 2 = plane is also 
not a glide plane, since following the reflection by the 



translation of (a x /2) (0, 1,0) to restore the z = a z /4 and 
z = —a z /4 layers would change the z — ±a z /2 lay- 
ers. The 6 rotations using x-, y-, and z-axes passing 
through the origin are readily seen to be symmetries of 
the layer, as are the reflections through the (110) and 
(110) planes. They take the point (x,y,z) into the 
points {{x, y, z) , {x, -y, -z) , {-x, y, -z) , (-x, -y, z) , 
(-y, x, -z) , (y, -x, -z) , (y, x, z) , (-y, -x, z)} which is 
the group D 2 d- This is a true point group, and the space 
group is therefore symmorphic. Spin-orbit effects come 
from terms linearly proportional to <7j, the spin operators 
that transform as pseudovectors, while the electric field 
E = (0,0, E z ) and in-plane momentum k = (k x ,k y ,Q) 
transform as vectors, the same as the coordinates. Our 
interest here is in combinations of these three quantities. 
Using the above list of operations, we find in zeroth order 
in E z that there is one invariant term in the Hamiltonian 
of the form k x a x — k y a y . In first order in E z , there is only 
the Rashba term E z (k x u y — k y a x ) (which is of course in- 
variant under all isometries). Either of these terms can 
be multiplied by any even function of E z , which is invari- 
ant under all the operations. Thus we find that when N 
is an odd number, j3 (N, E z , v) is an even function of E z 
and a (E z , N, v) is an odd function of E z . 



For even N, take the origin (0, 0, 0) at the center of 
a bond between atoms at positions (a x /8,a x /8,a z /8) 
and (—ax/8, —a x /8, —a z /8) . The origin is then a 
center of inversion. The rotation through 180° about 
the (HO) axis is a symmetry operation. The (110) 
axis is a screw axis since the 180° rotation about this 
axis must be accompanied by a translation through 
(a x /4, a x /4, 0) to be a symmetry operation. The same 
is true for the 180° rotation about the (001) axis. The 
8 operations of the factor group obtained by combin- 
ing these operations take the point (x,y,z) into the 
points {(x, y, z) , (-x, -y, -z) , {-y, -x, -z) , (y, x, z) , 
(y + a x /4, x + a x /4, -z) , (-y - a x /4, -x - a x /4, z) , 
(-x + a x /4, -y + a x /4, z) , (x - a x /4, y ~ a x /4, -z)}. 
Modulo translations, this is isomorphic to the group 
D 2 h- Because of the appearance of the translations, 
this is not a true point group and the space group 
is not symmorphic. Its action in the Hilbert space 
reduces in many cases to projective rather than faithful 
representations of Dih- However, this does not affect 
the symmetry analysis of the Hamiltonian. In zeroth 
order in E z , the group does not allow any combination 
of terms of the form ki<jj, since all such terms change 
sign under inversion. In first order, we again have the 
Rashba term E z (k x a y — k y a x ) . Again, multiplication 
of this expression by any even function of E z is permis- 
sible. Thus we find that when N is an even number, 
a (N, E Zl v) is an odd function of E z and (3 (E Z ,N, v) is 
an odd function of E z . 
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III. VALLEY AND SPIN ORBIT COUPLING 

We take into account the above symmetry arguments 
and employ a (5-functional approach for the interface- 
induced valley mixing 21 . In the envelope-function 
picture^, the Bloch periodic functions (valleys) at the 
interfaces \±)—u±k a e ±lkaL mix, resulting in valley split- 
ting, At, = A|<i>(zj)| 2 cos (koL — <fi v ), with A and <j) v be- 
ing phenomenological parameters, and k$ the conduc- 
tion band minima. |<E>(zi)| 2 is the zero-field value of 
the envelope function at the interfaces. Nestoklon et 
ali^ extend this work by introducing valley-orbit and 
spin-orbit mixing, as a spin-dependent reflection of the 
wavefunction at the interfaces, situated at z=z u ^ (see 
Fig. HJa)). Following their approach, we introduce 
next i?D and Hr as (5-functional perturbations in the 
lowest spinor-valleys functions at the interfaces. Wc 
consider the hight-symmetry directions x' || [110] and 



y' || [110], along which the spin eigenstates are paral- 
lel (see Fig. Q] (b)). In this rotated basis, we have: 
Hso = (a + ff)a x ik y i + {fi — a)a y >k x i. As mentioned in 
the text, a QW with N even is isomorphic with the D^k, 
which contains full inversion symmetry. Under the Dih 
operations, both Hr and Hry change sign: Ho — » —Ho] 
Hji — ► — Hr. However, for a QW with N odd, the 
point group is D2d, which transforms: Ho — > Ho', 
Hr — > —Hr. Taking into account these symmetry ar- 
guments, we introduce a (5-functional SOC, where the 
effective potential is at both interfaces situated in z u and 
Zd, as depicted in Fig. [TJa). The advantage of the 
5-functional is that it allows to consider the symmetry 
arguments in a simple way, but we shall recover the com- 
mon notation of |T]) and ([2|). The Hamiltonian in the 
(5-functional approach can be expressed in terms of the 
Pauli matrices in the valley-space, sf. 



H 



so 



/ 



Oy'kj 



'){aiS(z - z u ) - a*S(z - z d )) + (a X 'k y > + a y >k x >)(bi5(z - z u ) 



(-l) N b*6(z-z d ))], 
(3) 



where we have introduced and bi as parameters that 
determine the strength of Rashba and Dresselhaus SOC. 
We shall express those in terms of the commonly used a 
and (3. Eq. ([3]) can be divided in two contributions: the 
intra-valley SOC, which is Eq. ([3]) with i = 0, and the 
inter-valley SOC, for i — x,y: Hso = -ffsoi + Hso2- 

We first consider the intra-valley contribution, Hsoi (* 
= 0). so is then the identity matrix in the valley space, 
and thus, the spin-mixing parameters, ao and bo are real. 
Integrating ([3]) we recover the usual notation: 

Hsoi = [k x ,a y ,(p(E z ,N;0) + a(E z ,N;0)) 

+ k y ib x i (/3(E Z , N; 0) - a(E z , N; 0))] I, (4) 

where we have identified a and [3 in terms of the value of 
the envelope functions at the interfaces, <& u = Q(z u ) and 

a(E z ,N;0) = a = a |(|$„| 2 - |$d| 2 )|; 
P(E„N;0) ^/3 = 6oI(|^| 2 -(-in^| 2 )|. 

(5) 

Note that the parity is contained in the Dresselhaus 
term, which would be non-zero for N odd and E z — 
(or |$(z„)| 2 = |$(z d )| 2 )|). The Hamiltonian of Eq. Q 
mixes spins of the same valley, which we denote as 
{|+)} and {[—)}. Remember that in this basis, the two 
lowest eigenstates have even or odd parity 2 ^, ip e / — 
( e -i«V2. ijye^"/ 2 ), with rj = sgn{cos (k Na z /4)} deter- 
mining the parity of the ground state, and 4> v defined 
We include the spin degree of freedom to diago- 



,21 



{I t,+);| !,+>;[ T,-);l !,->}: 



-t- e -t(*<.-*i)/2 
_ e i{<j> v -<t>i)/2 
V T g*(<^+0l)/2 



with 4>i(ki 



arg{(l±i))sgn(/3 Ta )} 



I e -i(4>v+4>i)/2 \ 

-t- e -t(*<.-*i)/2 
^l)/2 
+0l)/2 

'(6) 

Note that 
=0, ±1, so 



/ 



{(p'\ d \a x ,\(p\ d ) = ±1, and {<p\ d \a y >\ip\ ti ) 
the eigen- vectors are eigenstates of a x < or to ay , depend- 
ing on the phase (f>% given by the direction of the in-plane 
wave vector, k. Along the x' or y' directions, the spin 
'up' and 'down' states split in both valleys by the same 
amount, [er-f — e jj = Ai oc 2|fc||(/3 ± «o)|- 



(a) 



Si 



\Na 7 




Dresselhaus 



FIG. 1: (a) Rashba and Dresselhaus spin- orbit fields for different 
directions ink space: along x' andy' directions, the eigenstates are 
parallel, (b) Schematic representation of a Si QW: the potential is 
considered in the interfaces at z u and z,j. 

We consider next the valley-mixing SOC, i?so2, 
which is Eq. ([JJ with i — x,y. The parameters, 



,b y =b=\b\e l - 



and a x =a v =a=\a\e 



i<f>a 



-b z =0) are 



nalize Eq. and find the eigen- vectors in the basis now complex, as they mix different valleys, and deter- 
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mine the valley-spin coupling strength. As previously 
noted by Nestoklon et ai, a Si QW possesses mirror ro- 
tation operation 54, resulting in a relative phase change 
of cf> a (4>p) for the Rashba (Dresselhaus) interaction at ei- 
ther interface. Combined with the absence (existence) of 
inversion center for even (odd) N , a change of sign is also 
observed in the Rashba (both Rashba and Dresselhaus) 
terms. Taking into account these symmetry arguments, 
the valley-mixing Hamiltonian has the shape: 



-ffso2 = / [(cr x >ky> - a y >k x <){a8{z - z u ) - a*S(z - z d )) 

+ {a x ,k y , + & y >k x ,)(bS(z - z u ) - (-l) N b*S(z - z d ))] dr. 

The valley-spin orbit Hamiltonian in the bi-spinor ba- 
sis for the lowest two valleys will thus consist of a 4x4 
matrix, which can be expressed in terms of the Pauli ma- 
trices for the valleys Sj and for the spin, <7j as: 



S02 



[a(E z ,N;i)(& X 'ky' - e y 'k x >) 



Cyf fax' )\ 



(7) 



with: 



a(E z ,N;x) = |a||(|$ u | 2 - |$ d | 2 ) cos0 Oq 

a(E z ,N;y) = |a||(|$ u | 2 + |$ d | 2 ) sin<fo Q 

/3(E z ,N;x) = |fe||(|$(z u )| 2 -(-l) A '|$(z d )| 2 )cos^ 

f3(E z ,N;y) = |fe||(|$(z„)| 2 + (- 1) ^ |<i>(^) | 2 ) sin^/3, 

where 0o = k$L and ^Oi = 4>Q ~ 4>i- We note that trans- 
lation of the vector r by a three-dimensional Bravais- 
lattice vector a results in multiplication of the Bloch 
functions |±) by the factors exp (±ikoL)2&, and thus a 
phase cf>o=koL appears in the valley-mixing terms. 



#S02 = y [<7y> kx>(j3z ~ Ot-z) + CT x >ky>(J} z + Otg)} , (8) 



with 



a|(|$ u | 2 cos0 Q - |^ d | 2 cos 0+) 

$„| 2 cos<^ - (-1)^1*003 0+), 



where we have defined: <j) i = 0oi ± <j>v Hence, we can 
write the SOC Hamiltonian in a compact way, merging 
Eq. (01) and Eq. ©: 



H so = Zy'k x ,J2 Si [P- n) (E z ,N)-a^\E z ,N) 

n i—0,z 

+ <W£ E § i [rt n) (E z ,N) + aM(Ez,N) 



n 2—0,2 



(9) 



Along k x i {k y i), the eigenvectors are eigenstates of ay 
(<j x i). The inter- valley term i = z, has a relative change 
of sign for the splittings of the spin 'up' and 'down' states 



in either valley, A2 oc ±2fc(/3 z ± a z ), so there is a valley- 
dependent spin splitting, as depicted in the inset of Fig. 
[21 |ef — £x| = |A 2 ± Ai|. From our numerical results, we 
observe that in general, IA2I > |Ai|, causing a reversed 
symmetry in the spin structure in the lowest two valleys. 
We also observe that higher order terms contribute to 
the valley-mixing SOC, as well as intra-valley SOC in 
the heterostructure case. Hence, we have generalized 
Eq. (|9|) by labelling the order n of the interaction, so 
far considered to zero order: aS n > oc (A:q^°^) 2 /|£„ — £o|. 
Note that the numerical results presented in this work 
correspond to a = J2 n and /3 = ^ n (3^- n \ 



IV. QUALITATIVE PICTURE 

Beyond the symmetry arguments, we can also analyze 
the functions a (N, E z , v) and j3 (E z , N, v) in qualitative 
terms. Let us first note that there are two distinct 
regimes for these functions considered in the E z — N 
plane. In the weak- field (WF), thick well regime, the 
wavefunction for the electrons in the lowest electric sub- 
band is spread throughout the well. The strong-field 
(SF), thin- well regime is reached when the wavefunction 
is confined near one interface and does not feel the other. 
In both the ideal and sandwich cases, the potential is 
rather flat in the interior of the well over a region whose 
extent ~ N and the confinement comes from relatively 
sharp interfaces. In this case, placing the classical turn- 
ing point in the middle of the well shows that the dividing 
line between the two regimes is described by 

(WF): E Z N 3 < 



n V 

= 1.5 x 10 11 — . 



(SF): E Z N 3 > 



miea z 
32ft 2 
mieo? z 



1.5 x 10 



11 V 



mi is the longitudinal mass. We need to consider the 
two sides of this line separately 

1. N dependence of a (E z , N, v) . For a the parity of 
N is not important. Let us define the lowest order term 
in E z for a (E Z ,N, v) as a 1 (N, v) E z (k x a y — k y a x ) . At 
first sight, the Rashba effect appears to be a bulk effect 
and therefore we expect a 1 (N, v) to be independent of 
N. However, the Ehrenfest theorem implies that the ex- 
pectation value of E z , which is proportional to the mean 
force, must vanish for any wavefunction bound in the 
z-direction. Thus in a continuum effective mass approx- 
imation the lowest-order term must vanish even though 
it is allowed by symmetry. Only when we put in inter- 
face effects and other atomic-scale effects will this term 
emerge. We shall assume that the extent of the inter- 
face in the z-direction is independent of N. If this is the 
case, then the probability to find the electron at the in- 
terface in the WF regime is ~ 1/iV, and we may expect 
a (E z , N, v) to be a decreasing function of N in the WF 
regime. In the SF regime (large TV) a becomes indepen- 
dent of TV for fixed E z since we only add layers that are 
unoccupied. 
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2. E z dependence of a (E z , TV, v) . We have seen that 
at small E z (WF) the dependence on E z is linear. At 
large E z for fixed TV (SF) the wavefunction is increasingly 
squeezed onto the interface and we may expect some con- 
tinued increase in a. Hence, we can set: 

a{E z ,N,v)=a 1 (N,v) ■ E z 

3. TV dependence of (3 (E Z ,N, v) . 

(a) TV odd. This is the only case for which 

(E z = 0, TV, v) ± 0. In the WF regime this field- 
independent term may be considered as a perturbation 
in 1/TV, since it is zero for even TV and the adding of an 
additional layer to make TV odd is the same as adding a 
term to the Hamiltonian whose matrix elements vanish 
as 1/TV. So we expect an initial decrease in the term as 
a function of TV. Again, a should approach a constant at 
large TV and fixed E z for the same reasons as in 1. 

(b) TV even. (E z = 0, TV, v) = 0. In the WF 
regime the field-independent term should converge to the 
result for even TV as TV increases, since they differ by 
terms of order 1/TV. The same holds for the SF regime. 

4. E z dependence of (3 (E Z ,N, v) . 

(a) TV odd. There is a constant term but no 
strong dependence on E z in the WF regime. In the SF 
regime the wavefunction is strongly confined to the inter- 
face. If we consider just a two-layer interface, there is a 
very strong orthorhombic anisotropy: the [110] and [110] 
directions are different, since the nearest-neighbor bond 
is in one of the two directions. Hence (3(E Z ,N, v) can 
be expected to be large, so (3 (E z , N,v) should increase 
strongly at large E z with fixed TV. 

(b) TV even. In the WF regime the symmetry is very 
important and (3 (E z , TV, v) is linear in E z . Again, in the 
SF regime we expect to converge to the odd TV result. 

Finally we consider the v dependence. The val- 
ley splitting A„ vanishes in the effective-mass con- 
tinuum approximation - it is due to interface effects. 
For E z = 0, the eigenfunctions must be even in z: 
■ip + (r) = F (z) <j) {f) cos k$z or odd in z: ip-{f) 

1 (z) <fi (r) sin fcoz, where F (z) is an even, slowly-varying 
envelope function and <f> (r) is an even (in z) Bloch func- 
tion, fco is the wavector of the conduction-band minima. 
For a well with smooth surfaces such as we consider here, 
E v has the order of magnitude ~ 1 — 10 meV and oscil- 
lates with thickness on the scale ATV = ir/2koa z and is 
proportional to 1/TV, as expected for an interface effect 
in the WF regime. In the SF regime, E z ^ and the 
eigenstates are no longer even or odd. A„ saturates for 
large TV at fixed E z , and its overall magnitude increases 
with E z , also as expected as the wavefunction is squeezed 
onto an interface. 

The oscillations with TV arise in the following way. 
Let V (z) be the confining potential and V (k) its 
Fourier transform: V (k) = J exp (ikz)V (z) dz. If 
we apply lowest-order degenerate perturbation the- 
ory for states in the two valleys, we find A v = 
\J d 3 r F 2 (z)cj) 2 (f)e 2lkoZ V (z)\ ~ \V(2k )\. As we 
change TV, V (z) has variation on the scale z = 4TVo, the 



separation between the two interfaces. |y(2fco)| then 
has constructive interference when 2ko x 4 (ATV) a = 2ir 
or ATV = n/Akoa. This ignores Umklapp, which will be 
present in the actual system, but is absent in the tight- 
binding approximation. It turns out that the dependence 
on the valley index can be quite dramatic. The valley 
states differ substantially right at the interface, where 
much of the spin-orbit effect arises. By the same to- 
ken, we can expect the same oscillations with TV that are 
seen in A„ to be present in Aso, the spin-orbit energy 
splitting. 

In actual heterostructures, the interfaces are not sharp. 
Ge is substituted for Si on randomly chosen sites, which 
will generally mean that the penetration length of the 
wavefunctions into the barriers varies randomly in x and 
y. All symmetries are violated by the disorder and it 
no longer makes clear sense to speak of even and odd TV. 
Furthermore, the free electrons come from dopants that 
create an electric field in the structure, so that the E z = 
limit is not accessible. It will probably be very difficult to 
observe the parity dependences that are predicted from 
symmetry arguments and also the oscillations. However, 
it may be possible to observe such effects in free-standing 
membranes. 



V. RESULTS FOR IDEAL CASE 

In this section we get the tight-binding results for the 
free-standing layer. Our approach to determining a and 
(3 will be to compute Ago along [110] and [110] using 
NEMO-3D for free-standing layers with varying thickness 
TV and applied electric field E z . To discriminate the 
Rashba and Dresselhaus contributions, we also compute 
the expectation value of a x i and a y > , which determines 
4>i of Eq. ([6]) and thus, the sign of ((3 — a). 
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FIG. 2: Splitting of levels due to spin mixing within the same 
valley and among different valleys. Note that A2 > Ai, in accor- 
dance with our numerical results. 

In Fig. [2] we show the four energy eigenvalues of the 
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lowest electric subband as a function of wavevector in 
the [110] direction for a 8 nm QW, which corresponds to 
N = 60, and E z — 0. The valley splitting is about 9 
meV, which is much larger than the spin-orbit splitting 
for all k. It is seen that Ago is linear in k at small 
k. Extracting A2 and Ai (see diagram in the inset of 
Fig. [2J we can then calculate [f3 + a) and (/3 — a) from 
the slopes of these lines, which allows us to determine 
a t [E z , N) and ft [E z , N) (see Eq. ©). 

We consider first an ideal Si QW, were the symmetry 
arguments can be directly applied. We impose hard 
wall boundary conditions at both interfaces and perform 
tight-binding calculations to obtain the eigenvalues and 
eigenvectors. As mentioned above, we expect the Rashba 
contribution to be linear in E z , a^(E z ,N) « ocqE z , and 
insensitive to N in the SF regime. The inset of Fig. [3] 
shows that this is the case: The intra-valley contribution 
of the Rashba SOC, aj, is plotted as a function of N. 
For large N , ao converges to a value and is insensitive 
to any changes of N, as corresponds to the SF regime. 
From this data, we find aj ~ 2 x 10 _5 nm 2 in the SF 
limit. The inset shows ao as a function of E z : for the 
regime of fields considered here, it appears to be linear. 
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FIG. 3: Rashba SOC contributions to the splittings of the first val- 
ley subbands in Si QWs as a function of thickness, N , for a unitary 
electric field, E z . Inset: ao as a function of E z for different QW 
thicknesses. 

For the Dresselhaus contribution, we consider first the 
N even case. We do not observe any SOC related terms 
at E z — 0, as expected: We recall that any term of the 
form fcjCTj is violated under symmetry operations. 
We expect the terms to be linear to lowest order of E z , 
which is indeed the case: (3q(z) — Po( z )Ez- The results for 
the intra-valley are shown in Fig. [3J The SF value for 
the intra-valley is 0q ~8x 10 -5 nm -2 . We note an abrupt 
change in /3q accompanied by a parity flip (depicted in the 
inset of Fig. , an event that has already been noted in 
literatur o 24 ! 25 ' 26 ! 27 . This reveals that (3q is more sensitive 
to higher energy contributions than ao. 

The inter- valley parameters are plotted next. We ob- 
serve again a linear behavior with electric field for both 
terms, as expected from our qualitative arguments. Fig. 
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FIG. 4: Intra-valley Dresselhaus SOC contributions to splittings 
of the first subbands in Si QWs as a function of N (N even). The 
insets show the valley-spin symmetry of the four lowest conduction 
subbands: a parity flip occurs at N ~ 20. 



|3] shows (a) a\ [a z = a\E z ), and (b) Pi {(3 Z = fi\E z ). 
Recall that under the D211 operations, both the Rashba 
(for any TV) and the Dresselhaus SOC (for even N) trans- 
form in a similar manner, hence a similar behavior is ex- 
pected with E z . We observe that both contributions 
exhibit oscillations as a function of N. These oscilla- 
tions are related to the valley-splitting oscillations, al- 
ready observed in literatur o 15 ' 21 i 23 ' 27 i 28 . For large N, 
the interaction converges to the SF limit value: a\ ~ 
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FIG. 5: Linear term of the inter-valley SOC contributions to split- 
tings of the first subbands in pure Si Q Ws, as a function of well 
width, (a) Rashba and (b) Dresselhaus (for even N) contributions. 
The lines are a guide to the eye. 



These results lead to two significant conclusions: (i) 
the Rashba inter- and intra-valley contributions are of 
the same order of magnitude, which would indicate that 
the bulk part of a predominates over the interface part 
in the SF regime. On the contrary, /3* is more than one 
order of magnitude larger than f3 z , indicating that it is a 
pure interface effect, also sensitive to higher order energy 
levels, (ii) We also find from Figs. E2[5]that the dividing 
line between the SF and WF regime is at E Z N 3 ~ 10 11 . 

Next, we consider the Dresselhaus contribution for N 
odd. Fig. 03 shows (3q as a function of E z for different 
N. We observe that for large N, the parity effect is less 
apparent, and /3q becomes independent of N, as predicted 
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from our qualitative arguments. In the WF regime, (3 Z 
is not as sensitive as (3$ to electric fields. We note that 
Po(E z — 0) is non-zero, as shown in the inset of Fig. [6] 
It presents strong oscillations with TV in the WF regime, 
due to mixing with higher energy states, and vanishes in 
the SF regime, as expected in the bulk limit for silicon. 
The overall Po(E z = 0) ~ 1/N dependence is also evident 
in this curve. 
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FIG. 6: Intra-valley Dresselhaus SOC contributions to splittings 
of the first subbands in Si QWs as a function of E z for odd N. 
The inset shows the zero-field value /3o(E z = 0) as a function of 
N. 



VI. RESULTS FOR HETEROSTRUCTURE 
CASE 

Having verified that the results are reasonable overall, 
we redo the calculations for a more realistic model of 
an actual heterostructure. First, we present results 
on a SiaGei-K/Si/Si^Gei-a; membrane with x = 0.5. 
The Si layer is surrounded by 28 layers of Si^Gei-x on 
both sides. This is sufficient to avoid surface effects 
and to confine the wavefunction in the Si QW for the 
electric fields presented here. Biaxial (shear) strain is 
minimized using VFF, which includes precise values of 
the elastic constants c,*, as defined elsewhere— The 
unit cell is chosen to have 24 atomic layers along the 
[010] and [100] directions, for which realizations of the 
substitutional disorder are averaged over, as noted in 
our numerical data. These calculations are much more 
time-consuming than those for the ideal layer, so fewer 
results are presented. We recall that calculations of 
this kind have not been done previously. Earlier work 
used the virtual-crystal approximation for the outer 
layers, which artificially preserves the symmetry. The 
reduction in symmetry can only increase the number 
of possible terms in the Hamiltonian, so not only the 
Rashba and Dresselhaus terms exist, but in principle all 
terms ki<jj can exist - strictly speaking, even k x and k y 
are no longer good quantum numbers. However, we 
shall take advantage of the approximate symmetry to 
present the results in the same way. 



The inter-valley Dresselhaus coupling constant for N 
odd is shown in Fig. [7J For large N, j3 z becomes lin- 
ear with E z , as it corresponds to the bulk limit. The 
behavior is quite similar to /3q. 
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FIG. 7: Inter-valley Dresselhaus SOC contributions to splittings 
of the first subbands in Si QWs as a function of E z , for N odd. 
The inset shows the zero-field value f3 z (E z = 0) as a function of 
N. 
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FIG. 8: Intra-valley contributions to SOC of a 

SixGei-x/Si/SixGei—x (x = 0.5) membrane as a function 
of E z , (a) Rashba and (b) Dresselhaus, for different number of 
atomic layers of Si, N . The inset of (b) sketches the structure, 
a SiGe/Si/SiGe suspended membrane. », N =20; I, N=36; A, 
N=72 ; o, N=108; A, N=1U; 

Fig. [8] shows the parameters a (a) and /3o (b) as a 
function of electric field, E z . We note that the dramatic 
dependence on the parity of iV is no longer present, as the 
disorder destroys their distinct symmetry properties. We 
also note that a is always non-zero, even for E z — 0. The 
intra-valley ao is non-linear in the WF regime, reaching a 
linear-in-i?2 value comparable to the one of Si QW in the 
SF regime (open triangles), ~ 1.9xl0 _5 nm 2 . This is 
consistent with the value obtained above for the Si QW. 

Fig. M (b) shows an overall 1/N dependency of Pi for 
E z = , consistent with the previous section results. 
For large N we observe that /?o is linear with electric 
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fields. We fit QW in the 10-30nm range and find 0k ~ 
38xl0 _5 nm 2 . We also find, for the SF regime, (3\ ~ 
58xl0" 5 nm 2 , with ft = (3}E Z . 

We have also studied the eigenvectors for each case and 
find frequent parity-flips (see inset of Fig. |4|) by varying 
TV or E z . As a consequence, both and a depend on TV 
even in the SF regime. 

Finally, we include the most experimentally relevant 
case: a Si^Gei-^/Si/Si^Gei-K grown on a Sii-^Gea; sub- 
strate, with x = 0.3. The in-plane lattice constant a\\ is 
now relaxed to the SiGe. We apply fixed boundary con- 
ditions for the lattice constant to the SiGe value at the 
bottom of the QW, and allow NEMO-3D to minimize 
the strain energy by varying the lattice constant along 
the z-axis. 




E z [xl0 6 V/m] E z [xl0 6 V/m] 

FIG. 9: Intra-valley contributions to SOC for a 

Si x Ge\- x /Si/Si x Ge\— x formed on a Si x Ge\- x substrate, 
with x =0.3. / The Rashba (a) and Dresselhaus (b) contributions 
are shown as a function of E z for different number of atomic 
layers of Si, N. The right vertical axis show the absolute energy 
values of the splitting, calculated for the Fermi level of a typical 
Si QW with n s = 4 X 10 11 cm -2 —. The inset of (b) sketches the 
structure: a Si layer is grown in SiGe substrate and capped with 
a SiGe layer. 

Fig. [5] shows the linear-in-fc SOC intra-valley con- 
tributions. The splittings are linear in E z in the SF 
limit, although parity flips and higher energy levels re- 
sults in small non-linearities. We note, however, that the 
zero field shows a visible splitting, even for wide QWs 
(TV ~140). We observe that ao is much smaller than 
in the previous samples, where silicon was the dominant 
component. For a typical Si/SiGe heterostructure with 
n s = 4 x 10 11 cm -2 — we have — 0.16nm _1 , at which 
the Dresselhaus-induced SOC splitting is ~ 1.25 fJ.eV, 
whereas the Rashba splitting is only 0.02 fieV. 

We show in table [T] the linear SOC coefficients for a 
20nm Si QW (TV ~140) in the three different structures 
considered in this work, for a typical electric field E z — 
10 7 V/m. a\ and ft 1 are in units of 10 _5 nm 2 . Note that 
the Dresselhaus is dominant in all cases. 



VII. CONCLUSIONS 

We have been able to extract the anisotropic-in-fc split- 
tings due to SOC using a sp 3 d 5 s* tight binding model ca- 
pable to take into account the interface effects to atomic 
scale. For a Si QW, we distinguish the TV even and odd 



TABLE I: eel and fij i= 0,z numerical results for 20nm Si QW 
in the three different fashions considered in this work, all in units 
of 10~ 5 nm 2 . 



a}(pl)[xW- 5 nm 2 ] 




aj 


Pi 


Pi 


'Pure' Si Membrane 


2.1 


1.3 


8.0 


154.5 


SiGe/Si/SiGe Membrane 


1.9 


2.5 


37.7 


58.5 


Si QW on SiGe Substrate 


0.7 


1.5 


30.6 


97.8 



cases, since symmetry operations over the Dresselhaus- 
type terms are fundamentally different: while for TV odd 
it appears at zero-order in E z , for TV even is linear in E z , 
to lowest order (as Rashba- type terms are). We have 
extracted the linear-in-TS z and linear-in-fc parameters for 
the TV even case. 

We also distinguish two regimes of operation for typ- 
ical wells: in the weak field regime, the splittings vary 
strongly as a function of TV. The intra-valley mixing 
components show roughly a 1/TV behavior, whereas the 
valley- mixing ones present also oscillations. On the con- 
trary, the splittings do not change with TV in the strong 
field limit. In the case of Dresselhaus, the dependency 
on TV is more prominent. Together with some oscil- 
lations observed for the intra-valley mixing, this would 
reveal that higher-energy states are also very important. 
We also observe a reverse spin structure in the spin-split 
valleys, a direct consequence of the inter-valley splitting 
being larger than the intra-valley. 

We find consistent results for Si QWs formed in SiGe 
heterostructures. We have also studied the lowest sub- 
band eigenstates, and found frequent parity flips by vary- 
ing TV or E z , suggesting a sample-dependent SOC. In 
accordance with ref.— , we find that the energy splittings 
due to Dresselhaus are in general larger than Rashba type 
ones: the Dresselhaus parameter (3 is almost one order of 
magnitude larger than the Rashba a for the ideal (pure 
Si QW) case, and roughly two orders of magnitude larger 
for the heterostructure (Si on SiGe substrate) case. New 
experimental data fits show indeed that the Dresselhaus 
term is larger than the Rashba one^i. We recall, how- 
ever, that throughout this paper the numerical data have 
been obtained for QWs with flat interfaces along a main 
crystallographic axis. More realistic samples include a 
small tilted angle with respect to a high symmetry axis. 
The numerical data are beyond the scope of this publi- 
cation, however, we expect the interface induced Dressel- 
haus parameter [3 to become smaller in tilted samples in 
a similar manner as valley splitting does^ 2 -. We expect, 
nevertheless, that a would remain of the same order of 
magnitude even in realistic tilted samples. Simulations 
carried out on a thin membrane confirm this hypothesis, 
although a more extensive study would quantitatively de- 
termine the SOC parameters. 
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